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1  INLEIDING 


Tot  de  insteUingen  noodzakelijk  voor  hci  maken  van  een  rijsimulatie  met  een  voertuig  op  de 
Universele  Dynamometerbank  (UDB)  bchoren  de  zogeheten  K-factoren.  Er  zijn  vier  K-factoren: 
Kj,  K2,  K3  en  K4.  Voor  de  bepaling  van  Kj,  K2  en  K3  is  op  het  Fysisch  en  Elektronisch 
Laboratorium  een  algoritme  ontwikkcld  (bet  zogeheten  schatteralgoritme)  dat  uit  een  op  de  weg 
gemeten  uitrolcurve  dczc  met  wrijvingscocfficienten  corresponderende  factoren  berekentCl.  K4  is 
bet  product  van  de  massa  van  bet  voertuig,  het  kwadraat  van  de  dynamische  wielradius  en  een 
schaalfactor.  Het  schatteralgoritme  kan  de  op  de  weg  gemeten  uitrolcurves  goed  benaderen:  de 
theoretische  uitrolcurves  die  bchoren  bij  de  berckende  K-factoren  vallen  praktisch  samen  met  de 
op  de  weg  gemeten  uitrolcurvcsl'l 

Het  doel  van  deze  evaluatic  was  na  te  gaan  of  voor  de  DAF  YA-4440  de  uitrolcurves  gemaakt  op 
de  UDB  overcenkomcn  met  de  wegmetingen. 
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2  THEORIE 


Aangezien  er  uitrolcurves  op  de  UDB  gemeten  zijn  met  66n  of  meer  K-factorcn  gelijk  gesteld  aan 
nul  kan  het  gebruikte  model  ten  behoeve  van  begripsvorming  vereenvoudigd  worden.  Omdat  de 
nadruk  bij  deze  evaluatie  ligt  op  het  rcproduceren  m.b.v.  K-factoren  van  een  uitrolcurve,  en  niet 
zozeer  op  de  waarde  van  fysische  cocfficienten,  wordt  de  fysische  bewegingsvergelijking 
herschreven  met  K-factorcn. 


2.1  De  bewegingsvergelijking 


Dc  bewegingsvergelijking  die  ten  grondslag  ligt  aan  de  Universele  Dynamometerbank  met 
fysische  codfficiSnien  f,.  {2  en  c^^  luidt: 

M  =  m  r2  (1)'  -I-  m  g  r  fj  -t-  m  g  r2  f2  CO  +  A  p  r3  c»  (1)2  /  2  (1) 

waarbij; 

M  het  op  de  aandrijfas  aangrijpcnde  moment  voorstelt; 
m  de  massa  van  het  vocrtuig  voorstelt; 
r  dc  dynamische  wiclradius  voorstelt; 

(0  de  hoeksnclhcid  van  cen  wicl  van  het  voertuig  voorstelt; 

(O'  de  eerste  afgcleidc  van  dc  hocksnelheid  van  een  wiel  naar  de  tijd  voorstelt; 
p  de  luchtdichtheid  voorstelt; 

g  de  graviialicversnclling  voorstelt; 

A  het  cffectieve  oppcrvlak  van  het  voertuig  voorstelt; 
f,  de  snclhcidsonafhankclijkc  wrijvingscocfficient  voorstelt; 

f2  de  lincair  snclhcidsafhankclijkc  wrijvingscoCfficicnt  voorstelt; 

c,^  de  luchtwccrstandscocfficiem  voorstelt. 


TNO-rapport 


Pagina 

7 


De  conversieformules  om  K,,  K2  en  K3  uil  fj,  {2  en  le  berekenen  volgen  hieronder: 
mgr 

K,  =  -  f, 

12 

50  m  g  r2 

K2=  -  f2 

3 

500  A  p  r3 

K3=  - c, 

3 

Voor  K4  geldt; 

K4  =  2*mr2/5 

Hieruit  volgt  dc  bewegingsvcrgclijking  met  K-factoren: 

M  =  (5/2)  K4  (O'  +  1 2  K,  +  (3/50)  Kj  (0  +  (3/1000)  K3  0)2  (2) 

Tijdcns  uilrollcn  is  M  allijd  gelijk  aan  nul; 

5  K4  (O'  =  -  24  *  K,  -  (3  /  25)  *  K2  (0  -  (3  /  500)  K3  0)2  (3) 


2.2  De  opiossingen  van  de  differentiaalvergelijking  voor  bijzondere  gevallen 

Voor  ijkingsdoclcinden  worden  in  dit  hoofdstuk  dc  opiossingen  van  de  differentiaalvergelijking 
gegeven  voor  verschillcndc  combinatics  van  K|,  K2,  K3  en  K4.  cOq  staat  in  de  opiossingen  steeds 
voor  de  bcginhocksnelheid  ten  tijde  van  dc  start-  of  referentietijd  t  =  0. 
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De  bewegingsvergelijking  voor  (Kl,  K2,  K3,  K4)  =  (^,  0,  0,  >0)  is; 

5  K4  0)'  =  -  24  *  K,  (4) 

De  oplossing  van  deze  vergelijking  is: 

cod)  =  -(24/5)  *  (K,/K4)  *  l  +  (Oo  (5) 

Merk  op  dat  co  lineair  afhankelijk  is  van  dc  lijd. 

Indien  alleen  K3  nul  is  luidi  dc  bewegingsvergelijking; 

5K4d(0/dt  =  -24*K, -(3/25)*K2(0  (6) 

Dc  oplossing  is: 

cod)  =  ((Ci+C2*(Oo)/C2)*EXP(-(C2/C4)t)  -(C,/C2)  (7) 

waarbij: 

cOq  siaat  voor  dc  bcginsnclhcid; 
c,  =  24  K, 

C2  =  (3/25)  K2 
C4  =  5  K4 


Merk  op  dat  co  exponcnticcl  afhankelijk  is  van  t. 
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Voor  (Ki,  K2,  K3,  K4)  =  (?i0, 0,  >0,  >0)  luidl  de  bewegingsvergelijking: 

5K4(O’  =  -24*Ki-(3/500)K3{o2  (8) 

De  oplossing  van  de  bewegingsvergelijking  is  nu  afhankelijk  van  hei  teken  van  Kj.  Omdat  een 
negatieve  Ki  een  constante  acceleratietenm  voorstclt,  wordt  alleen  de  oplossing  voor  positieve  Kj 
gegeven: 

0)(t)  =  sqrt(Ci/C3)*tan(arctan(sqrt((c3/c,)*a)o))  -  sqrt(c,*C3)*t/C4)  (9) 

waarbij; 

c,  =  -24  K, 

C3  =  3  K3/5OO 
C4  =  5  K4 


Merk  op  dat  co  niei-lineair  afhankelijk  is  van  i. 


Sndhdd  (km/u)  SneHieid  [km/u] 
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3  CONSISTENTIE 


Om  de  UDB  te  controleren  op  consistentie  zijn  er  voor  een  aantal  instellingen  van  de  K-factoren 
uitrolcurves  gemaakt  op  de  UDB  met  en  zonder  DAF  YA-4440. 


Fig.  3.1.  Drie  UDB-uitrolcurves  met  (K,J<2J(3.K4)=(36.20,170,754)  met  DAF  YA-4440. 


Tl|d  (si 


Fig.  3.2.  Drie  UDB-uitrolcurves  met  (KiJ<2,KjJ<4)=(23,95, 1 13,754)  zonder  voertuig. 


TNO-rapport 


Pagina 

11 


Deze  en  andere  hier  niet  getoonde  grafieken  geven  aan  dat  de  UDB  met  en  zonder  DAF  YA-4440 
consisiente  uitrolcurves  produceeri;  in  de  grafieken  vallen  immers  de  uiirolcurves  prakiisch 
samen.  De  uitroltijden  verschillen  maximaal  3%.  Figuur  3.2  loom  het  grootste  geconstateerde 
verschil  in  de  uitroltijden. 
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4  IJKING 


Wanneer  theorctische  uitrolcurves  vergeleken  worden  met  UDB-uitrolcurves  gemaakt  zonder 
voertuig  voor  (Ki,K2.K3,K4)  =  (>0,0,^,>0).  dan  blijkt  dat  er  een  verschil  tussen  de  uitrollijden 
bestaat.  Daartoe  een  tweeial  voorbcclden: 


Fig.  4.1.  Theoreiische  en  UDB-uitrolcurve  zonder  voertuig  met  (K]J(2J(s,K4)=(30,0, 0,754). 
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Fig.  4.2.  Theoretische  en  UDB-uitrolcurve  zander  voertuig  met  (K],K2J<^3d(4)=(30, 0,1 50,754). 

Dc  UDB-uitrolcurve  blijki  systematisch  cen  langere  uitroltijd  te  hebben  dan  volgens  de  theorie 
zou  moctcn.  Dc  verschillcn  in  dc  uitroltijd  blijkcn  bij  uitrolcurves,  waarvoor 
(Ki,K2.K3,K4)  =  (>0,0,0,>0),  proportioned  te  zijn  met  de  uitroltijd.  Voor  een  aantal  verschillende 
waarden  voor  Kj  zijn  de  verschillcn  in  de  uitroltijden  berckend  tussen  UDB-uitrolcurves  gemaakt 
zondcr  voertuig  en  theoretische  uitrolcurves.  De  resultaten  (K4  is  constant  754)  zijn  in  label  1 
weergegeven. 

Omdat  K,  klcincr  is  dan  K4  heeft  een  klcine  absolute  verandering  in  Kj  grotere  gevolgen  voor  de 
UDB-uitrolcurve  dan  een  kleinc  verandering  in  K4.  Het  is  dus  voor  de  hand  liggend  om  een  offset 
te  bcpalcn  in  Kj  en  niet  in  K4. 

Uit  de  bcwegingsvergelijking  (4)  voor  (K,,K2,K3,K4)  =  (>0,0,0,>0)  volgt: 


AK,  =  ((OE,nd-cOo)*K4*(-5/24)/At 


(10) 


waarbij; 
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cOq  de  startsnelheid  voorstelt; 

^Eind  de  eindsnelheid  (snelheid  waarop  metingen  gestaakt  wordt)  voorstelt; 

At  het  verschil  in  de  uitroltijdcn  om  le  decelereren  van  (Oq  tot  a>Eind  voorstelt; 
AKi  is  het  verschil  tussen  de  ingestelde  Kj-waarde  en  de  weiicelijke  Kj-waarde. 


Tabel  1 


K, 

ingesield 

Uitroltijd 

UDB-bank  zondcr 

voerluig 

[s] 

Uitroltijd 

theorie 

[s] 

Verschil 

uitroltijdcn 

At 

[s] 

AK, 

30 

188.5 

208.8 

-20.3 

-2.9 

50 

113.1 

120.6 

-7.5 

-3.1 

75 

75.4 

78.7 

-3.3 

-3.1 

100 

56.5 

58.3 

-1.8 

-3.0 

150 

37.7 

38.5 

-0.8 

-3.1 

200 

28.2 

28.8 

-0.6 

-3.6 

Mcrk  op  dat  de  onbctrouwbaarhcid  van  AKj  toeneemt  naarmate  At  afheemt,  terwijl  de  absolute 
fout  in  At  in  dc  orde  van  0.2  -  0.3  s  blijft.  Dc  enigzins  afwijkcndc  waarde  voor  AKj  voor  Kj  =  200 
t.o.v,  de  gemiddclde  waarde  van  AK,  is  dan  ook  niet  verontrustend. 


Uit  tabel  1  kan  geconcludcerd  wordcn  dat  K,  een  offsetwaarde  heeft  ter  grootte  van  3.0. 

De  theoretische  uitrolcurves  met  (K|,K2,K3,K4)  =  (ki-3>0,0,k3>0,k4>0)  vallen  samen  met  de  op 
de  UDB  gemeten  uitrolcurves  zondcr  vocrtuig  met  (Ki,K2,K3,K4)  =  (k,>0,0,k3^,k4>0)  (zie  fig. 
4.3  en  4.4).  Eventuele  afwijkingen  in  K2  cn  K3  blijken  in  dit  specifieke  geval  nauwelijks  van 
bclang  te  zijn. 


Sndheid  (km/til  Smlhek)  (km/u] 
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—  Zonder  voertuiQ 


Fig.  4.3.  Theoretische  uitrolcurve  met  (KjJ(2J(3,K4)=(27, 0,0,754)  en  UDB-uitrolcurve  met 
(K ,J<2,Ki,K4)=(30,0,0,754)  zonder  voertuig. 


0.00  20  00  40  00  60.00  80  00  100  00  120.00  140.00 

Tl|d  Is) 

Fig.  4.4.  Theoretische  uitrolcurve  met  (Kj,K2J(3,K4)=(27,0, 150,754)  en  UDB-uitrolcurve  met 
(Kj,K2J<s,f(4)=(30.0,l50,754)  zonder  voertuig. 
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5  WEG-  EN  GESIMULEERDE  UITROLCURVES 

5.1  Vergelijking  theoretische  wegcurve  en  gesimuleerde  uitrolcurve 


Om  een  uitrolcurve  gemaakt  op  de  weg  te  simuleren  met  de  UDB  werden  de  met  het 
Schattcrprogramma  voor  de  DAF  YA-4440  berckende  K-factoren  ingesield,  waarbij  K2  gelijk 
gesteld  werd  aan  nul  en  K)  vermeerderd  werd  met  3.  Hierbij  wordt  verondersteld  dat  de 
snelheidsafhankelijke  wrijvingsterm  (K2)  vooral  een  gevolg  is  van  interne  frictie  in  de 
aandrijftrein  van  het  vocrtuig  zelf  en  dus  niet  gesimuleerd  dient  te  worden  op  de  bank,  indien 
datzelfde  vocrtuig  daarop  getest  wordt. 

In  figuur  5.1.1  wordt  de  optimalc  curve  gcbaseerd  op  UitrollAl*)  getoond.  De  door  het 
schatteralgoritme  bepaalde  K-factor  in.stellingen  zijn  (K,.K2,K3,K4)  =  (22,94,114,754).  Eveneens 
in  de  figuur  is  de  berckende  UDB-uitroIcurve  gemaakt  met  de  gecorrigecrde  K-factoreninstelling 
(K,,K2,K3,K4)  =  (25,0,1 14,754). 


90  00 

80  00 

70.00 

60  00 

^  50  00 

I  40.00 
cS 

30  00 
20.00 
10.00 
000 


\ 


- - ypg  YAJ1440 

Weg  beiekend  YA-4440 


0.00  20.00  40  00  6000  80  00  lOOOO  120.00  140.00 


Tiid  [sj 

Fig.  5.1.1 .  Berckende  UitrollA  en  de  UDB-uiirolcurve  met  (KjJ(2,K^J(4)=(25,0,l  14,754)  en 

DAF  YA-4440. 
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In  figuur  5.1.2  worden  de  bcrckende 
(K„K2,K3,K4)  =  (23,96,1 11,754)  en  de 

(K,,K2,K3,K4)  =  (26,0,1 1 1,754)  getoond. 


curve  voor  de  wegmeting  Uiirol2Al') 
bijbehorende  UDB-uitrolcurve  met 


80.00 

70.00 

60.00 

—  50  00 

1 

40.00 

*«) 

30  00 

20.00 

10.00 


UDB  YA-4440 

Wefl  berekend  YA-4440 


0.00  - - - , - 1 - - - 1 - 1 

0.00  20.00  40  00  60  00  80.00  100.00  120.00  140.00 

TIJd  [s] 

Fig.  5.1.2.  Berekende  Uitrol2A  en  de  UDB-uitrolcurve  met  (KiJ(2,KjJ(4)=(26,0,1 1 1 ,754)  en 

DAF  YA-4440. 


De  gesimuleerde  uitrolcurves  bcnadercn  de  berekende  uitrolcurves  goed. 


5.2  Invioed  kleine  verandering  K| 


Om  de  invioed  van  een  kleine  verandering  in  K,  na  te  gaan  werden  verschillende 
UDB-uitrolcurvcs  gemaakt  met  voor  K,  afwijkcndc  waarden. 
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60  8 

Tlldlsl 


Fig.  5.2. L  Berekende  UitrollA  en  UDB-uitrolcurves  met  (KjJ(2J(^,K4)=(25-i-A,0J  14,754)  met 

A=-1.0.1  enDAF  YA-4440. 


Uitrol2A 

Delta=-1 

Della^O 


De 


Fig.  5.2.2.  Berekende  Uitrol2A  en  UDB-uitrolcurves  met  (KiJ<2J(j,/(4)=(26+A,0,l  11,754)  met 

A=-l.0.1  en  DAF  YA-4440. 


Hieruit  blijki  dat  een  klcinc  vcrandcring  in  K|  gecn  grote  gevolgen  heeft  voor  de  iigging  van  dc 
gcsimuleerde  uilrolcurve  t.o.v.  dc  uil  dc  wcgmcUngen  berekende  curve. 
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53  Gesimuleerde  en  gemeten  uitrolcurves 

Tot  slot  worden  gesimuleerde  en  op  de  weg  gemeten  uitrolcurves  met  elkaar  vergeleken.  Enige 
voorzichtigheid  is  op  zijn  plaats  aangezien  de  K-factoren  afgeronde  waarden  zijn  en  in  de  praktijk 
gebleken  is  dat  verschillende  sets  K-factoren  eenzelfde  uitrolcurve  goed  kunnen  beschrijven. 

90  00  n 

-  UDB  YA-4440 

1  >  Weg  gemeten  YA-4440 
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Fig.  5.3.1.  Gemeten  UitrolIA  en  de  UDB-uitrolcurve  met  14,754)  met 

DAF  YA-4440. 
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6  CONCLUSIES 


De  instellingen  voor  de  K-factoren  voor  de  DAF  YA-4440  kunnen  bcpaald  worden  m.b.v.  hel 
schatteralgoritme,  en  geven  een  gocdc  overeenstemming  tussen  wegmetingen  en 
simulatormetingen.  Voor  K]  is  hicrbij  een  offsetwaarde  ter  grooite  van  3  in  acht  genomen.  K2  is 
hierbij  op  de  UDB  gelijk  gesteld  aan  nul.  Onzeker  is  of  K2  gelijk  aan  nul  stellen  voor  alle  typen 
voertuigen  geldl. 

De  correclies  op  de  resultaien  van  hel  schatteralgoritme  voldoen  voor  de  DAF  YA-4440.  Er  kan 
echter  nict  gcgarandeerd  worden  dat  de  coirecties  in  het  algemeen  voldoen  voor  hel  uitvocren  van 
rijsimulaties. 

Verdere  verbetering  is  waarschijnlijk  mogelijk  indien  de  K-factoren  ingesteld  kunnen  worden  met 
groiere  relatieve  nauwkeurigheid. 

Naar  aanlciding  van  deze  resultaien  wordt  een  andere  kaart  in  de  UDB  geinstalleerd,  waarmee  de 
K-factoren  met  hogere  relatieve  nauwkeurigheid  ingesteld  kunnen  worden. 
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